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In this paper we present fits to charmless hadron ic B decay data from the BaBar, Belle and 
i-^3 ' Cleo experiments using models by Beneke et. al. anc j Ciuchini et When we include 

the data from pseudo-scalar vector decays (PV) the current experimental results favour the 
O 1 inclusions of a so-called "charming penguin" term. We also present fit results for the Unitary 

Triangle parameters and the CP violating asymmetries 
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•i-h , 1 Introduction 

X: 

A wealth of experimental data on hadronic charmless B decays has become available from the 
BaBar and Belle experiments. In this paper we present an analysis, based upon QCD factori- 
sation, of data on hadronic charmless B decays. We also investigate the potential contribution 
to the decay amplitudes of b quark annihilation and of so-called charming penguins. The data 
we attempt to fit includes pseudo-scalar (tt and K) and vector (p, u, K* and 4>) mesons. This 
is an extension of an earlier study^. 

The decay amplitude can be described by: 

A(B -> MM) = -t J2 E V? <M 1 M 2 \Q l \B> fact (1) 



p=u,ci=l..lo 



with X p the appropriate CKM matrix factor for the process. The coefficients Oj can be written 
as: m = dij + ctiji. 

The dij coefficients are similar to the ones of naive factorisation with vertex and penguin 
corrections included. We assume that all the light mesons have the same spatial wave function. 
As a consequence the an coefficients are universal. The diji coefficients contain the hard 
spectator interactions. They contain an end-point divergence from the low energy contributions 
to the loop integrals, which is parametrised by the complex parameter Xh- The magnitude of 
Xh is expected to be between and 3. 



B meson decay can also be initiated by the b quark annihilating with its partner. Contri- 
butions to the decay amplitude should be small, of order Aqcd/Mb, compared to the leading 
b quark decay. Because of the heavy b quark mass it is expected that perturbative QCD calcu- 
lations will give a reliable estimate of the contributions of annihilation to the decay amplitude. 
Apart from the electroweak meson decay constants, Jb, fir, f P etc, there is a low energy con- 
tribution coming from the loop integrals. This is parametrised by the complex parameter Xa- 
Formula for the annihilation contributions are given by BBNS , we use the extension to pseudo- 
scalar vector decays by Du et al.^. 

Charming penguins orginate from diagrams with a charm quark loop. According to Ciuchini 
et al. they can be enhanced^. They behave similar to the 04 coefficients in signs and Clebsch- 
Gordan coefficients and take the following form: 

A^ q = ^ [-v: q v ub {h - a gim ) - v: q v qb Pi] (2) 

2 Fitting Method 

We have attempted to fit the theoretical expressions for branching ratios with the available data 
as averaged by the Heavy Flavour Averaging Group^. Measured branching ratios for seventeen 
channels are shown in Table ^ We take the measured branching ratios to be the mean of the 
B and B decays. The CP asymmetries are not included in the fit. Their measurements are 
not always consistent between the experiments and the errors are large. Therefore we prefer to 
compare the measured results with the prediction from the fit. 

For convenience we assign to each channel (/11/12) a number a. The statistical and systematic 
errors have been combined into a single error o a . The systematic errors are small in general and 
we ignore any correlation between the channels. We then write for a x 2 function 

X 2 (Pi) = ^2 [(Bra(-Pi) — Br Q (exp)) /cr a } 2 + additional constraints. (3) 

a 

Br a (Pj) are the theoretical branching ratios in terms often parameters Pi, i = 1, . . . , 10 which we 
take to be the three Wolfenstein CKM parameters {^4, p, rj} and the seven soft QCD parameters 
{R^, , Fk, Fk, A p , A^, Ak*}- For the fit to the charming penguin model we introduce four 
more terms: P 1 and P GIM and their phases. In this case we fix A to the world average on 0.82 
and keep RZ' fixed at 1.0. The well known decay parameters {/„-, fx, f P , fu>, ftp, Ik*} are held 
at their mean values and the Wolfenstein CKM A parameter is taken to be A = 0.2205. For the 
divergence parameters we take \Xh\ = \Xa\ = 2, which leads to reasonable values of e.g. sin 2/3. 
Additional terms were included in the x 2 to take into account experimental and theoretical 
constraints from outside the data on B decay branching ratios. We search for a minimum of \ 2 
as a function of the Pi. We use the MINUIT^ program to minimise the x 2 ■ The theoretical 
branching ratios and contributions of the individual channels to x 2 based on these best fit values 
are given in Tabic ^ 

The values of the best fit parameters are shown in Table [2] together with our estimates of 
the errors. These errors are of course highly correlated. A plot of the error matrix ellipse for 
the Wolfenstein parameters p and r/ is shown in Fig. ^ For both models, the results in Tabled 
for the best fit values of the various form factors lie within the spread of theoretical estimates. 
The x 2 /dof is 34/16 for the BBNS fit and 16/12 for the charming penguins. 

3 Conclusions 

The first conclusion would be that the factorisation approach works well. Most of the branching 
fractions in Table Q are predicted correctly by both models. The fitted parameters in Table [21 



Table 1: Measured branching ratio Br(exp), experimental error a (errors added in quadrature), theoretical branch- 
ing ratio for best fit parameters and contribution to x 2 f° r various B decay channels. Units are 10 -6 . 
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look reasonable for both fits. The x 2 f° r the charming penguin fit is significantly better. This is 
due to the poor fit of BBNS for the decay modes with a K* meson. Also the experimental value 
for ujK° is not easily accommodated within the BBNS model. Both models have a problem 
fitting the 7r°K° mode. Figure H shows the position of the apex of the Unitarity Triangle for 
both fits. The results agree with each other on the angle 7, but give a larger value than the fits 
from other measurements. The angle (3 agrees well, although it has to be pointed out that for 
the BBNS the value of (3 is sensitive to the choice for Xa- 

Regarding the asymmetries it is too early to come to a conclusion. In many cases the 
experiments do not agree, in others the errors are so large that a meaningful discrimination is 
not possible. The asymmetries are in principle very sensitive to the different models and with 
improved statistics could become the final test of factorisation. 



Table 2: Best fit values and one-standard deviation errors for both methods. F Wi k and A PtUlt K * are the transition 
form factors, F B ^ M (0). R^' K are the chiral enhancement factors, which are nominally power supressed, but in 

practice O(l). 
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Table 3: CP asymmetries from BaBar and Belle and fit predictions. Only statistical errors are included. 
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Figure 1: Result for the Unitarity Triangle fit for BBNS (left) and charming penguins (right). The shaded area 
shows the 3<r allowed region for the apex of the Unitarity triangle. The data point shows the fit result from the 
Unitarity Triangle fit from other measurements taken from C KM fitter^ . 
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